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1. INTRODUCTION 

Induction motors (or asynchronous motors) now account for about half of total electricity use in 
industrialized countries. Furthermore, A.C. motors consume a significant amount of power in the agricultural 
and industrial sectors. The quantity of energy consumed by a motor over its life cycle is calculated to be 
60-100 times the motor's initial cost [1]. Technological characteristics enjoyed by this type of induction 
motors squirrel cage, such as low maintenance cost, high reliability, high performance, and a good starting 
torque. Also, the type of rotor has the easy and most rugged construction imaginable and is nearly 
indestructible. Therefore, squirrel cage induction motor is one of the most versatile electric machines of 
previously and today, where being form in approximately 70% of global industry and also consuming around 
50% of the energy generated on the planet and other applications about machines [2], [3]. 

The induction motor has consisted of a laminated stator and rotor core. The classic design induction 
motor stator has one set of 3-phase windings that are spatially distributed uniformly around the stator and is 
driven by a sinusoidal balanced polyphase excitation that drives currents through the stator windings, in 
which symmetrical three-phase voltages displaced by 120 (electrically) are supplied [2]. While the stator of a 
6-ph of induction motor (I.M) comprises two groups for different 3-ph windings that are offset by a particular 
angle from one another. It has axial rotor bars that are often cast in aluminum or copper [4]. Electric current 
passes through the stator conductors, producing a magnetic field that rotates at a synchronous speed. The 
single cage rotor is the most common rotor bar type. Single layer windings have been used in small (AC) 
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motors, and a single-layer winding is an example of a three-phase concentric winding [5]. Mathematical 
models have been created and analyzed by two concept three-phase and six-phase induction motor [6]. The 
criteria for selecting the starting points for each phase were proposed as a general formula for determining 
the number of slots available for an N-phase alternating current machine design [7]. In industrial drives, six- 
phase induction motors have numerous benefits over traditional 3-phase motors including reduced magnetic 
flux harmonics, increased stability, reduces torque pulsations, and lower power levels for the static converter 
[8]. It is important to establish models for these phenomena to gain a deeper understanding of them. The 
numerical electromagnetic field measurement is good method for designing and calculating induction motors 
or electrical machines. Other researchers proposed a method to describe a 3-ph linear I.M to determine the 
different values employed in its /phase duplicate circuit by a DSP system [9]. It has found that the air gap in 
linear induction motors (LIMs) is much larger than in rotary induction motors (RIMs). The secondary 
parameters have determined using two methods based on MATLAB-Simulink [10]. 

For analyses of electrical machines, the finite element method (FEM) uses the subdivision of a large 
problem domain into more simple sections called finite elements and variational methods to solve the 
problem by reducing a related error function [11]-[13]. MATLAB simulation based on the symmetrical 
component theory, the estimated performance of a three-phase IM at steady-state, such as under/over and 
unbalanced/balanced voltage disturbances at different values of VUF, has been analyzed [14], [15]. The 
T-slip characteristics and S.state power losses have also been measured and plotted at various degrees of 
voltage unbalance. The modern programs help design the electrical machines on the reducing cost of the 
application and which lead to achieving the best results for the design before building real the structure of the 
machines [14]. 

COMSOL multiphysics is one of the most successful simulation programs used in design the 
induction motor and it's used in other applications [16]. It has been used in this current work to analyze 
induction motor parameters. The "rotating machinery" interface in COMSOL multiphysics 5.5 has used to 
run a frequency domain and non-linear simulation of an induction machine's dynamic activity action. The 
two-dimensional FEM model is linked to electrical circuits by connecting the physics interfaces "rotating 
machinery" and "electrical circuit" [17]. Also, the COMSOL multiphysics has included the multi-time 
analysis in the same program like time-dependent, frequency-domain [18]. The two-dimensional model in the 
COMSOL program is a perfect method for studying induction motor behavior under different situations, 
speed, torque, and current values are depending on the precision of identification and calculation of induction 
motor parameters, geometry, and material properties [19]. 

FE program, traditional analysis calculation, and actual measurement have been used to look into 
the flux density in the air gap, torque output of the dual 3-ph motor [18]. Trapezoidal wave phase current has 
used by investigating the air-gap field density of a 6-phase induction motor [20]. The induction motor has 
represented as a two-dimensional model with specified out-of-plane thickness. As a quasi-3D model. The 
COMSOL AC/DC module's rotating machinery, the magnetic interface is used for stationery and time- 
domain modeling [21], [22]. The dependent variables in this physical interface, which satisfies Maxwell's 
equations, are the magnetic vector potential and magnetic scalar potential [21]. The magnetic interface of the 
nodes ampere's law has used in the rotating machinery to solve the magnetic field in the motor [23]-[25]. The 
main problem of the 3-phase induction motor is the beginning of the start operation torque and distribution of 
magnetic flux density. Therefore, this present paper suggested enhancing the performance design of a 
3-phase induction motor squirrel cage by COMSOL multiphysics. 

The contribution of this work to the first time has used the flux density results and change of the 
core material in the same design by COMSOL software program. In this paper has been changed the 
magnetic flux density (MFD) for the three-phase and six phases of the squirrel cage induction motor and 
compares the two cases above. In addition, two materials have used frequency-domain analysis in the design 
of a six-phase induction motor. The first material is aluminum (AL), and the second cast iron (CI) material. 
These materials have applied to the rotor part of the motor on six phases of induction motor to obtain the best 
results. FEM has been used to analyze the result of torque (N/m) with omega speed (rad/s) to know which 
material is better than the design of the six phases in the induction motor. The values of the parameter of the 
electrical conductivity, relative permeability, and relative permittivity have taken from the default of the 
same program (add materials) of the COMSOL. 


2. RESEARCH METHOD 

A 3-phase induction motor's performance can be improved by modifying the MFD and the core 
materials of the rotor part. This design requires high accuracy for the same parameters of the motor to keep 
the comparison between three and six phases of design. 
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2.1. Induction motor construction 
The induction motor is made up of two parts: 
a) The stator, which is the outer component 
b) The rotating component 
The stator and rotor are both constructed of: 
a) The electric circuit, which is often built of insulated aluminum or copper winding, carries current 
b) The magnetic circuit, which is typically constructed from a series of laminated silicon steel laminations 
pressed together to produce a cylindrical magnetic circuit to carry MFD current 
There are many important equations to consider when designing an induction motor stator [26]-[31]. 
For full pitch, pitch factor (K, = 1). The motor output power (Pout) is calculated as in (1): 


Pout = HP X 746 (1) 


where (HP) stands for horse power in a motor. The input apparent power (KVA;) of a motor can be 
determined as in (2): 


nx pf (2) 


where (n), (pf) are a motor's efficiency and power factor respectively. The synchronous speed of the motor 
in revolutions per second (N,, ) can be calculated as in (3): 


fe 
s = 37 3) 

where (f) is the frequency of supply. The output coefficient (C,) is calculated as (4): 
Co = 1.117? Bay ac 107 (4) 


where Bay represents the average flux density and is the ampere conductor per meter. The stator's diameter 
and length are as in (5)-(8). 


psiz, = AV (5) 
TOOK Neg 
DSi?L x P 
DSi? = ——__ 6 
: mx L/t (6) 
DSi = VDSi3 (7) 
DSi? L 
= (8) 
DSi? 
where: 
DSi = the stator's inner diameter, 
L = length, and 
T = pole pitch. 
The magnetic flux (Ø) can be determined as (9). 
pe e laa (9) 
P 
The number of turns per phase (T;,;,) is calculated as (10): 
V, 
T ph (10) 


ph 4A4XK, Xf xO 
where (Vpn) is the voltage per phase. The total number of conductors (Z+) is calculated as in (11): 
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where (n) denotes the total number of phases. The number of conductors per stator slot (ZS,) is as (12). 


aes (12) 


J= ph (13) 


where (Acu) is the area of one conductor. And the rated per phase current (Ipa) is calculated as in (14). 


_ KVA; (14) 
ph nx Von 
The area of conductors per stator slot (AS,) equals: 
AS, = ZS; X Acu (15) 
and the area of the stator slot (AS,) is calculated as in (16): 
Hie (16) 
S 1-—spf 
where: (spf) denotes the space factor. The diameter of one stator slot (DS,) is given by (17). 
AS 
DS, =2x z (17) 
T 
The stator's outer diameter (DS,) is calculated as in (18). 
DS, = DS; + 2 x (2 x DS;) (18) 


The air gap length (Lg,) is calculated as in (19), 


Lg, =0.2+2x./DS; XL (19) 


and the rotor's outer diameter (DR,) is determined as in (20). 
DR, = DS; — 2 x Lg, (20) 


2.2. Induction motor with three phases 
The following equations would apply to 3-phase induction motors. The quantity of slits/pole every 
phase () is determined as in (21): 


Sse = 3x (21) 


where (P) is the number of stator slots and is the number of poles. The slot electrical angle (a) is calculated 
as in (22). 


_ 180 xP 


ae =— (22) 
S 


Furthermore, the distribution factor (K4) equals: 
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sin = x Se) 
Kg = ——— E (23) 
S X sin- 


spp 
The winding factor (Kẹ) was calculated as in (24). 

Ky = Kp X Ka (24) 
The total number of conductors () is calculated as in (25). 


Z =2Xx3 X Tpn (25) 


The rated per phase current (lpn) is determined as in (26). 


KVA; 
lpn = 
3 x Von (26) 
2.3. Induction motor with six phases 
A 6-phase induction motor has a different stator than a 3-ph induction motor. The following are the 


equations for a 6-ph induction motor. The quantity of slits/poles per phase (Sspp) determined as (27). 


Ss 
Ssp = Exp 


(27) 
While the distribution factor is one for a six phase, four pole, twenty-four slot stator, the (Sspp) is one. 
Ka = 1, and Ky = 1 
The total number of conductors (Z+) is calculated as in (28). 
Zp =2X6XThynr (28) 


The rated per phase current (lpn) is determined as in (29). 


i i. 29 
ene eae 29) 


The main dimensions of model geometry are two-dimensional and is depicted in Figure 1. 


Figure 1. Depicts the dimensions of an induction motor 


Int J Elec & Comp Eng, Vol. 12, No. 1, February 2022: 62-72 


Int J Elec & Comp Eng ISSN: 2088-8708 o 67 


2.4. Result of stator design 

The following are the requirements for the induction motor to be designed: 2 horse power (HP), 
4 pole (P), supply frequency of 50 hertz (f), supply voltage 400 volt (VL), efficiency of 80% (n), full pitch, 
0.825 full load power factor (pf), unity ratio of pole length to pole pitch (L/t), 0.25 space factor (spf), 
24 stator slots (Ss), average flux density (Bav) of 0.44 Wb/m2, electrical loading (ac) of 18000 ac/m. 

A 20 SWG enameled copper coil, Dc (conductor diameter) = 0.914 mm, and Acu = 0.656118 mm? 
conductor area are used for three phase winding. Table | displays the outcome of designing the stator of a 
3-phase IM. A 22 SWG enameled copper coil, diameter of the conductor (Dc) = 0.711 mm, and area of the 
conductor (Acu) = 0.397035265 mm? are used for the six-phase winding. Table 2 shows the outcome of 
designing the stator of a 6-phase IM. 


Table 1. Design details for a two-horsepower 3-phase induction motor stator 


Parameter Value Unit 
DSo 157 Millimeters 
DSi 112.6879193  Millimeters 
DSs 10.83174626 Millimeters 
DSc 8.123809697 Millimeters 
DRo 112.0879193 Millimeters 
Tph 316 Turns 
Iph 3.262903794 Ampere 

J 4.9730405 A/mm2 
ZSs 79 Conductor 
Lgl 0.3 Millimeters 


Table 2. Design details for a two-horsepower stator for a 6-phase induction motor 


Parameter Value Unit 
DSo 159 Millimeters 
DSi 111.3931845 = Millimeters 
DSs 11.80250482 Millimeters 
DSc 8.851878614 Millimeters 
DRo 110.7931845 Millimeters 
Tph 310 Turns 
Iph 1.631451897 Ampere 

J 4.109085618 A/mm2 
ZSs 155 Conductor 
Lgl 0.3 Millimeters 


The per-phase current in a six-phase stator is half that of a three-phase stator; with a lower current. 
A smaller cross-section area for wire has used, so the normal wire gauge has been modified from 20 to 22, 
and the current density in the six-phase windings is around 82 percent of the three-phase windings, a ratio of 
current density as in (30). 


Jsix—pnase 


Current density ratio = = 0.83 (30) 


three—phase 


The number of turns and conductors measured as an integer number. After rounding, the number of 
turns for each phase of the 3-phase motor is 316 turns, while for each phase of the six-phase motor is 
310 turns; this small difference is primarily due to the winding factor. As shown in the preceding tables, the 
number of conductors per slot varies greatly, with 79 conductors per slot for three-phase stators and 
155 conductors per slot for six-phase stators. However, since the cross-section area of 22 SWG wire is less 
than that of 20 SWG wire, the number of conductors per slot ratio is as in (31). 


LSS ix hase 
Conductor per slot ratio = ——~"“"_ = 1.96 (31) 
Sthree—phase 


3. SIMULATION RESULTS 

The design of induction motor has been simulation by COMSOL multiphysics 5.5. The design has 
taken to compare between the three and six phases to know how the performance-enhancing by two 
simulation results, first the distribution of flux density to it, and the second changing of the core materials in 
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the FD to know the relation between torque and speed. Figure 2 shows the distribution of MFD in the 3ph 
induction motor design. Figure 4 shows the distribution of MFD of three-phase induction motor when it has 
taken the effect of rotor part rotational considered. Also, it observed the weak flux density on the distribution 
of rotor part, and the directional of rotation is anti-clockwise. Figure 5 shows the distribution of MFD of six- 
phase induction motor when it has taken the effect of rotor part rotational considered. Also, it observed the 
regularity flux density on the distribution of rotor part, more uniform intensity, and the directional of rotation 
is anti-clockwise. 


Figure 2. Three-phase MFD distribution in induction Figure 3. Six-phase MFD distribution in induction 
motor motor 


Omega(61)=1200 rad/s Surface: Magnetic flux density norm (T) B 


Contour: Magnetic vector potential, z component (Wb/m) 
mm 


-0.4 
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Figure 4. Distribution of MFD in the rotor part for three-phase induction motor 
when FEM has taken considered 
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Figure 5. Distribution of MFD in the rotor part for six-phase induction motor 
when FEM has taken considered 


Figure 6 shows the mesh configuration established of the six-phase induction motor in the air gap 
way with free triangles, and the color blue is shown slots of conductor winding in the stator part. Figure 7 
shows the relation between torque (N*m) and angular velocity Omega (rad/s) when has been chosen the 
frequency domain analysis in a 3-ph induction motor. Also, it detected the starting torque beginning and 
ending in the negative value. Because the field strength is low and, therefore the torque is small and starts 
with a negative value. It is worth mentioning the AL core material is employed in the rotor part of the 
3-phase induction motor because it does not show different results when the materials are changing. 

Figure 8 shows the relation between torque (N*m) and angular velocity Omega (rad/s) when has 
been chosen the frequency domain analysis in a six-phase induction motor. Also, it observed the starting 
torque beginning is positive because of the increase in the flux density and ending in the negative value at the 
steady-state. The material has used in the rotor part is aluminum material. 
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Figure 6. Mesh configuration of the six-phase induction motor 
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Figure 8. The 6-phase induction motor's torque and angular velocity relationship 
when the aluminum (AL) material 


Figure 9 shows the relation between torque and angular velocity Omega when has been chosen the 
frequency domain analysis in a six-phase induction motor. Also, it observed the starting torque beginning is 
positive because of the increase in the flux density and ending in the negative value at the steady-state. The 
material has used in the rotor part is CI material. It is observed the starting up to 0.8 and steady-state 
approximately to equal -1. 
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Figure 9. The relation between torque and angular velocity omega of the six-phase induction motor when the 
cast iron (CI) material 
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4. CONCLUSION 

The scope of the paper is computer simulation help for the design of induction motors in their basic 
form useful to design adjustment. For the same design has a centered on two situations, first the MFD 
simulation three-phase and six-phase and the changing of the core materials of the rotor part for the induction 
motor. Second, changing of materials has been used with frequency domain in COMSOL program for the 
same design. The results showed the superior performance of the first situation for the six-phase MFD for 
time dependence and frequency domain when taken the FEM in considered. Also, the results have been 
shown in the second situation the positive effect of CI compared with AL on the performance in the relation 
between torque and angular velocity omega. 
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